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ABSTRACT 

We report here on multifrequency radio observations of the pulsed emission from PSR 
B1259-63 around the time of the closest approach to its B2e companion star. There 
was a general increase in the pulsar's dispersion measure and scatter broadening, and 
a decrease in the flux density towards periastron although changes in these parameters 
were seen on timescales as short as minutes. The pulsed emission disappeared 16 days 
prior to periastron and remained undetectable until 16 days after periastron. 

The observations are used to determine the parameters of the wind from the Be 
star. We show that a simple model, in which the wind density varies with radius as 
r~^, provides a good fit to the data. The wind is clumpy with size scales < 10^" 
cm, densities of ~10^ cm~^ and a velocity of ^--^2000 km s~^ at a distance of 20 - 50 
stellar radii. We find a correlation between dispersion measure variations and the pulse 
scattering times, suggesting that the same electrons are responsible for both effects. 
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1 INTRODUCTION 



The pulsar PSR B1259— 63 is a member of a unique binary 
system. Discovered using the P arkes telescope in a su rvey of 
the Galactic plane at 1.5 GHz ( Johnston et al. 1992a| ), it was 
shown by Johnston et al. (1992b) to be in a highly eccentric 
3.4-yr orbit with a lOth-magnitude Be star, SS 2883. 

Extensive observations of this system have been made, 
particularly near the previous periastron passage on 1994 
January 9. These include observations at radio wavelengths 
(llohnston et al. 199^; iMcClure-Griffiths et al. 1998[ 



et aL 1998 
1994) , in the X- 



19991), opticaT 
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ray (ICominsky, Roberts fc Johnston 



Kas£^et^L 1995 \_ Hirayama et al. 1996 ) and in 7-rays ( Grove 



et al. 1995; Hirayama et al. 199£). These observations have 



Johnston ct al . 



led to a better understanding of th e physical conditions in 



the wind and disc of the Be star (Pavani, Arons fc Kasp: 



199£) 



1995; Pavani fc Arons 1997; Johnston et al. 199 



1993 ; king fc Cominsky 1994|; iMelatos, Johnston fc Mclrosq 



Ball ct al 



The understanding of the system is thus as follows: SS 
2883 is a (main-sequence) B2e star with mass ~10 Mq and 
radius ~6 R0 . Its characteristic emission disc extends to at 
least 20 R, in the optical, similar to the distance between 



the pulsar and Be star at periastron. It is likely that the disc 
is highly tilted with respect to the orbital plane; the orbital 
plane is itself at 36° with respect to the line of sight. Due 
to the misalignment of the emission disc with the orbital 
plane, the pulsar goes behind the disc (with respect to the 
observer) near T -16 and emerges from the front of the disc 
near T +14. (Here, T is the epoch of periastron). The disc 
material is extremely dense near the stellar surface but falls 
off rapidly as a function of distance. The polar wind, by 
contrast is ~10* times less dense than the disc material at 
the stellar surface and has an density profile. 



The radio observations of the pulsed emission from the 
1994 periastron were rather crude in the sense that fre- 
quency agility was limited, with frequency switching requir- 
ing manual receiver changes and hence only possible on a 
~1 day timescale. Furthermore, only relatively few observa- 
tions were obtained before and after the eclipse of the pulsar. 
In 1994 there was no real-time pulse folding software which 
limited our best observing strategy in that the off-line data 
reduction significantly lagged behind the observing. How- 
ever, these observations showed dispersion measure (DM) 
and scintillation parameter changes on timescales of at least 
as short as one day and demonstrated that these changes 
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were h ighly asy mmetric with respect to periastron (John- 



al. 1996) 



ston et 

In this paper we report on observations of the pulsed 
radio emission made around the latest periastron passage of 
the pulsar on 1997 May 29. Section 2 describes the observa- 
tions. In Section 3 we present the results and in Section 4 
discuss their implications. 



2 OBSERVATIONS 

In 1996, a major upgrade was made to the focus cabin of 
the Parkes radio-telescope which allowed several receivers 
to be placed on a movable platform. This enables the choice 
of receiver to be made from the control room and allows for 
frequency changes in only a few minutes. 

We observed PSR B 1259-63 regularly with the 64-m 
Parkes radio tele scope at frequencies of 1.4, 4.8 and 8.4 GHz 
(VVex et al. 1998). The observations were intensified around 
the periastron passage of the pulsar on 1997 May 29. At 
each frequency the receiver consists of a dual-channel, cryo- 
genically cooled system sensitive to two orthogonal linear 
polarizations. 

Data were recorded simultaneously using two recording 
systems. The Caltech Fast Pulsar Timin g Machine (F PTM) 
is a correlator capable of fast sampling (Navarro 1994). The 



signals from the two polarizations are down-converted to 
intermediate frequency and 2-bit digitised by the FPTM 
which computes correlation functions (CFs) with 128 lags 
over a 128 MHz bandwidth. The system can be used either 
to record full Stokes parameters for one 128 MHz bandwidth 
or total intensity for two independent bandwidths. For the 
higher frequencies we always recorded the data in full Stokes 
mode, with the band centered at 4800.25 and 8400.25 MHz. 
In the 20cm band, we either recorded the full Stokes parame- 
ters at 1351.0 MHz or two independent bandwidths centered 
at 1274.75 and 1396.0 MHz. The CFs are then hardware in- 
tegrated at the apparent pulsar period (~47 ms) and the 
data recorded every 60 s. A graphical output displays the 
pulse profile every integration cycle. Thus, changes in ob- 
serving strategy could be made in real-time. 

Simultaneously, we also recorded data using a filterbank 
detection system in an identical fashion to that described 
elsewhere (e.g. Johnston et al. 1996). The signals from the 
receivers were down-converted to an intermediate frequency 
and then passed into a filterbank which consisted of 64 chan- 
nels, each 5 MHz wide for a total of 320 MHz. Full Stokes 
parameters were recorded every 0.6 ms and written to mag- 
netic tape for off-line analysis. The centre frequencies were 
1374, 4748 and 8348 MHz. 

For fiux density calibration, observations were made of 
Hydra A at 1.4, 4.8 and 8.4 GHz once every few days at the 
telescope. These measurements yielded system flux densities 
of ~28, ~100 and ~100 Jy at the three frequencies. Before 
every observation of the pulsar, a pulsed noise source, in- 
jected at 45 degrees to the signal probes, was recorded. This 
allows for calibration of the system gains and phases. 



2.1 Derivation of the Dispersion Measure 

Over a period of a few hours we obtained observations at 
5 different frequencies: 1.2, 1.4, 1.5, 4.8 and 8.4 GHz. The 



three frequencies in the 20-cm band were obtained simul- 
taneously. The observation would then be followed by one 
at either 4.8 and 8.4 GHz before again observing the 20- 
cm band. Typical observations were 1-hr in duration. Pulse 
time-of-arrivals were obtained in the standard manner (e.g. 
Wex et al. 1998). The DM was obtained by assuming the 
DM contribution from the 1.5 kpc path through the inter- 
stellar medium is 147.6 cm~^pc and does not change with 
time. Then, the extra DM contribution from the wind of the 
Be star could be computed by the relative differences in the 
residuals between the different frequencies. Having closely 
spaced frequencies in the 20-cm band allows any ambigui- 
ties to be resolved (i.e. where the ADM could impose more 
than 1 phase turn on the residual). Typical errors using this 
method are small (~0.1 cm~''pc), much smaller than those 
obtained simply from fitting the delays in frequency-phase 
space over a single observing band. 



2.2 Derivation of the Scattering Time 

Scatter-broadening of a pulse can be approximated by con- 
volving the intrinsic pulse shape with a truncated exponen- 
tial function. The situation for PSR B1259-63 is somewhat 
complicated by th e presence of two strong puls es, only 20 
ms apart in time (Manchester & Johnston 1995). For large 



scattering times, the scatter tail of the first pulse overlaps 
the second pulse and this must be taken into account in the 
derivation of the scattering time. We used the Levenberg- 
Marquart method (LMM) to perform a least-squares fit to 
the observed profile given a 'standard' (unscattered) profile. 
The fit returns the 4 free parameters : amplitude, baseline 
and decay time of the exponential and the offset in time be- 
tween the profiles. Generally, for scattering times less than 
7 ms, the fits are good. Above this, however, the pulses tend 
to be weak, and the extremely scattered and weak pulses 
make the fitting less robust. This is reflected in much larger 
errors bars on the fltted parameters. 



3 RESULTS 

The broad overview of the observational data is as follows. 
A change in the DM over and above that from the ISM 
was flrst observed at T -48. Over a period of ~20 days, 
both the DM and the pulse scattering showed an overall 
increase, and strong hourly variations. The pulsed emis- 
sion was last detected on T -17. Observations made on 
T -15 showed no pulsed emission at any frequency, a re- 
sult confirmed with independent ob servations on the Aus - 
tralian Telescope Compact Array (Johnston et al. 1999). 



The pulses were re-detected at Parkes on T -1-16, following a 
non-detection two days earlier. Post-periastron, no change in 
DM or pulse-scattering could be measured, although signif- 
icant changes in rotation measure and polarization proper- 
ties were recorded. Significant fiux variability was seen both 
prior to and following the eclipse. The general trend though 
was towards lower fluxes in the few days before and after 
the eclipse. 

Figure 1 shows the DM changes and the scattering time, 
r, as a function of epoch prior to periastron. Figure 2 shows 
a shorter time span of only 8 days prior to the eclipse. 

The figures show a wealth of detail. Points to note are 



© 0000 RAS, MNRAS 000, 000-000 



The 1997 periastron passage of PSR B1259—63 3 



: I I f 
I 1 ^ I 1- 



Days from Periastron 



Figure 1. Dispersion Measure changes and derived scattering 
times as a function of time prior to periastron. The errors on the 
DM are of order 0.1 cm~^pc and are too small to be seen on the 
plot. 




Time from 1997 May 8 (hours) 

Figure 3. Top panel : Flux density of PSR B1259-63 on 1997 
May 8. The observing frequency is 1.4 GHz and the fluctuations 
are broadband. Lower panel : Scattering time in ms for those 
integrations with largo enough s/n to obtain reliable estimates. 
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Figure 2. As Fig. 1, showing a period of only 8 days just prior 
to the eclipse of the pulsar. 



that the DM changes on short timescales (~1 hr or less) by 
~1 cm~"^pc. There is a slow increase in the DM from T -50 
until -27. The DM then slowly drops, before rising steeply 
again at T -17. The pulsar was undetected following these 
observations until T +16. Both the daily variations and the 
drop in DM towards T -20 are evidence for the clumpy 
nature of the wind. 

Flux density variations in PSR B 1259-63 occur on a 
number of timescales. The diffractive scintillation timescale 
at GHz frequencies close to periastron is expected to be very 
short (<C 1 second) and hence unobservable. The refractive 
timescale is probably of the order of minutes. One also ex- 
pects changes in the flux density due to free-free absorption 
as the pulsar moves behind clumps in the Be star wind. In 
Figure ^, we show flux variations as a function of time start- 
ing on 1997 May 8.3 (T -21). Each point is an integration 
over 2 minutes of data. These variations are broad band (at 
least across the 300 MHz observing bandwidth) at 1.4 GHz. 



At higher frequencies, the pulsar is too weak to allow reliable 
flux densities to be obtained in a few minutes. 

This extreme flickering is only observed on T -21. Ob- 
servations made the previous day and on the following three 
days show only slow changes in flux density over hours and 
a much lower modulation index, consistent with the refrac- 
tive scintillation timescale. T -21 corresponds exactly to the 
first en hancement in the flu x density seen in the unpulsed 
source ( Johnston et al. 1999| ) . This early increase in the flux 
density subsequently disappeared and does not fl t well with 
the model of the unpulsed source (Ball et al. 1999), although 
Johnston et al. (1999) have argued that it is likely to be the 
result of the first interaction between the pulsar wind and 
the stellar disc coupled with optical depth effects. 



4 DISCUSSION 

At T -20, the pulsar is ~50 stellar radii from the Be star. 
Assuming a wind density at the stellar surface of 7x lO^cm""^ 
which varies with radius as r~^, the line-of-sight density to 
the pulsar should be ~3xl0^cm~''. However, one cannot 
expect this wind to be homogeneous and there are likely to 
be significant density variations because of turbulence. 

We can use the equation derived by McClure-Griffiths 
et al. (1998) to relate the electron density in the scattering 
material, 5n^, to the scattering parameters. 



Snt = 4.0 Ai/T 



5/6 



n-5/6 

kpc 



11/3 



,2/3 
tOjpc 



{1 + x) 



5/6 



Here, is the scintillation bandwidth, D the distance to 
the pulsar and v the observing frequency, lo is the outer scale 
length, AL is the thickness of the scattering screen and x 
is the ratio of the observer- screen distance and the pulsar- 
screen separation. We simplify this by using D = 1.5 kpc, 
1/ = 1.4 GHz and assuming lo ~ AL, j: S> 1 to obtain 



5ni = 9.8 Ai/, 



-5/6 



5/6 



(2) 
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Figure 4. Scattering time versus change in DM at a given epoch. 
The straight hne shows t <x DM^^/^. 



Converting the scintillation bandwidth, /S.v to the scattering 
time, T, in milliseconds via 2iitAi' = 1 yields 

5nl = 1.43 X 10^ T^"" x^"" cm"^ 



(3) 



of ~ 8.4 X 10 



12, 



cm, 



5.5 X 10*), a scattering time of 3 ms gives 5n1 of 7 x 
cm~® and hence rie of ~ 8.3 x lO^cm"^ assuming a fully 



Assuming a distance to the screen, r 
(i.e X 
10" 

modulated wind. This is remarkably close to the expected 
value obtained from a simple scaling of the electron density 
in the wind and shows that the scattering is caused by the 
global properties of the fully modulated wind. The size of 
the scattering disc, given by {2TcrsY^^ is rather small, only 
3.9 X 10^° cm. 

These values are consistent with the DM measurements. 
A mean wind density of ~ 8.3 x lO^cm"^ along a pathlength 
of roughly 2rs yields a ADM of 4.5 cm~'^pc, close to the 
observed values. We also note that to double the observed 
scattering tail requires only a change in of ~30 per cent 
and that this change affects the DM at the ~1 cm~"^pc level. 
Such changes are seen on roughly a 2-hr timescale 

From equation 3 above it is evident that t (x n, 
whereas DM oc rig. Providing the same electrons are respon- 
sible for the scattering and the dispersion we should observe 
that r oc DM^^^^ Fi gure 2 shows the expected correlation 
between dispersion measure and scattering time in that they 
both rise and fall together. In Figure ^ we show a plot of 
scattering time versus ADM. The scatter in the diagram 
leads to a large uncertainty in any formal straight line fit. 
For illustrative purposes, we show a line with slope 12/5 as 
derived above, which is in satisfactory agreement with the 
data. 

Consider the free-free optical depth along the line of 
sight, given by 



12/5 



Tft = 8.2 X 10"^ T~ 



dl 



(4) 



where T is the temperature in Kelvin, u the observing fre- 
quency in GHz, d the distance through the absorbing mate- 
rial in pc and Ue the free electron density within the clump 



For an e lectron temperature in the wind of ~ 10 K (Water: 



IQSeQhe wind parameters described above yield ts ~ 0.6 
at 1.4 GHz, fully consistent with the lower flux densities 



recorded at these epochs. Flux density variations are seen 
both on timescales of a few minutes and on a roughly similar 
timescale to the DM and scattering variations. The longer 
timescale is probably due to variations in the free-free opti- 
cal depth, either because of temperature variations or small 
variations in the electron density. The short timescale vari- 
ations generally have a low modulation index and are likely 
to be refractive scintillation. Given a timescale of ~200 s for 
the variations and knowledge of the size of the scattering 
disc (3.9 X lO^^cm), the dominant velocity in the system is 
then ~2000 kms~^. Such a velocity is consistent with UV 
observations of Be star winds in general. 

In summary, we have shown that a plausible model for 
the wind of the Be star can account for the dispersion mea- 
sure and scattering data before periastron and is also con- 
sistent with the variability timescale and the flux density 
variations. We stress that we are discussing the spherical, 
high- velocity, low-density wind from the star. The pulsar en- 
counters the circumstellar disc material at T -18, the density 
of the disc is so large that the pulsed flux becomes eclipsed 
shortly thereafter through a combination of free-free absorp- 
tion and extreme pulse broadening and remains so until the 
pulsar re-emerges from the disc at T +16. Post-periastron, 
the pulsar lies between the observer and the Be star. The 
path length through the stellar wind is much less than pre- 
periastron and the DM and scattering variations are hence 
much smaller. 
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